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Abstract 

In this article we consider the formation and evolution of black holes, especially 
those in binary stars where radiation from the matter falling on them can be seen. 
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00 . We consider a number of effects introduced by some of us, which are not traditionally 

^^ I included in binary evolution of massive stars. These are (i) hypercritical accretion, 

which allows neutron stars to accrete enough matter to collapse to a black hole 

0> . during their spiral-in into another star, (ii) the strong mass loss of helium stars, 

5^ '. which causes their evolution to differ from that of the helium core of a massive star. 

(iii) The direct formation of low- mass black holes (M ~ 2M0) from single stars, a 

consequence of a significant strange-matter content of the nuclear-matter equation 



O \ of state at high density. We discuss these processes here, and then review how they 

-*^ ' affect various populations of binaries with black holes and neutron stars. 

^ ' We have found that hypercritical accretion changes the standard scenario for the 

evolution of binary neutron stars: it now usually gives a black-hole, neutron-star 
^ , (BH-NS) binary, because the first-born neutron star collapses to a low-mass black 

H ' hole in the course of the evolution. A less probable double helium star scenario 

has to be introduced in order to form neutron-star binaries. The result is that low- 
mass black-hole, neutron star (LBH-NS) binaries dominate the rate of detectable 
gravity-wave events, say, by LIGO, by a factor ~ 20 over the binary neutron stars. 

The formation of high-mass black holes is suppressed somewhat due to the infiu- 
ence of mass loss on the cores of massive stars, raising the minimum mass for a star 
to form a massive BH to perhaps 80 Mq. Still, inclusion of high-mass black-hole, 
neutron-star (HBH-NS) binaries increases the predicted LIGO detection rate by 
another ~ 30 %; lowering of the mass loss rates of Wolf-Rayet stars may lower the 
HBH mass limit, and thereby further increase the merger rate. 

We predict that ~ 33 mergers per year will be observed with LIGO once the 
advanced detectors planned to begin in 2004 are in place. 

Black holes are also considered as progenitors for gamma ray bursters (GRB). 
Due to their rapid spin, potentially high magnetic fields, and relatively clean en- 
vironment, mergers of black-hole, neutron-star binaries may be especially suitable. 
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Combined with their ten times greater formation rate than binary neutron stars 
this makes them attractive candidates for GRB progenitors, although the strong 
concentration of GRBs towards host galaxies may favor massive star progenitors or 
helium-star, black-hole mergers. 

We also consider binaries with a low-mass companion, and study the evolution 
of the very large number of black-hole transients, consisting of a black hole of mass 
~ 7 Mq accompanied by a K or M main-sequence star (except for two cases with a 
somewhat more massive subgiant donor). We show that common envelope evolution 
must take place in the supergiant stage of the massive progenitor of the black hole, 
giving an explanation of why the donor masses are so small. We predict that there 
are about 22 times more binaries than observed, in which the main-sequence star, 
somewhat more massive than a K- or M-star, sits quietly inside its Roche Lobe, 
and will only become an X-ray source when the companion evolves off the main 
sequence. 

We briefly discuss the evolution of low-mass X-ray binaries into millisecond pul- 
sars. We point out that in the usual scenario for forming millisecond pulsars with 
He white-dwarf companions, the long period of stable mass transfer will usually 
lead to the collapse of the neutron star into a black hole. We then discuss Van den 
Heuvel's "Hercules X-1 scenario" for forming low-mass X-ray binaries, commenting 
on the differences in accretion onto the compact object by radiative or semiconvec- 
tive donors, rather than the deeply convective donors used in the earlier part of our 
review. 

In an appendix we describe the evolution of CygX-3, finding the compact object 
to be a black hole of ~ 3 Mq, together with an ^ 10 Mq He star. In a second 
appendix we do the accounting for gravitational mergers and in a third one we 
show low-mass black-hole, neutron-star binaries to be good progenitors for gamma 
ray bursters. 

Key words: black hole physics — stars: binaries: close — - gamma rays: bursts — 
stars: neutron — gravity waves — accretion 
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1 Introduction 



The fate of massive stellar cores, both in single and binary stars, has many 
observable consequences, both for what types of compact object may be found 
in what type of binary, and for the formation rates of all types of compact- 
object binary. We have discussed various aspects of this problem in previous 
works, and here give an overview of all these together, applying the same set of 
principles to all and obtaining a consistent picture of the evolution of massive 
stars and binaries. 

The best-known compact-object (i.e., neutron star or black hole) binaries are 
the binary neutron stars. They are key testing grounds of general relativ- 
ity, and the usually favored gravity-wave source for LIGO. Until recently the 
theoretical formation rate of binary neutron stars gave at least one order-of- 
magnitude higher rate than was arrived at empirically by extrapolation from 
observed binary neutron stars. Because there are few binary neutron stars, 
and even fewer dominate the empirical estimates, the latter are frequently re- 
vised. The recent doubling of the estimated distance to PSR 1534-1-12 [2] has 
lowered the empirical birth rate significantly, widening the gap. 

A solution to this discrepancy comes from combining the strange-matter equa- 
tion of state, which results in a relatively low maximum mass for neutron stars, 
with hypercritical accretion [1] . In the standard scenario the first neutron star 
formed spirals into the other star, in a phase of common-envelope evolution. 
Bethe & Brown [1] argued that when a neutron star spirals into a red giant 
with a deeply convective envelope, it accretes matter at a very high rate of 
up to 1 Mq yr~^. Photons are trapped in the flow and carried adiabatically 
inwards to the surface of the neutron star [3] . The latter is heated to T ~ 1 
MeV, temperatures at which neutrino emission can carry off the thermal en- 
ergy. Hence, the Eddington limit of MEdd ~ 1-5 x 10^^ Mq yr^^ does not apply. 
As a result, the neutron star accretes about a solar mass of material and col- 
lapses to a low-mass black hole. Only if the two stars are initially so close in 
mass that at the time of spiral-in the first supernova has not yet exploded (i.e. 
the object that spirals in is still a helium star) a binary neutron star is formed. 
The sum total of binary neutron stars and black-hole, neutron-star binaries 
is almost the same as what was found for binary neutron stars in previous 
estimates, but now the binary neutron stars are only a small fraction of the 
total. The result is that an order of magnitude more black-hole, neutron-star 
binaries than binary neutron stars are formed. Together with the fact that 
the black holes are somewhat more massive than neutron stars, this implies 
that binaries with black holes should play an important part in mergers pro- 
ducing gravitation waves. They may also be good candidates for producing 
gamma-ray bursts. 



No low-mass black-hole, neutron-star binaries have been observed. This is due 
to the fact that the one neutron star in them is unrecycled, hence is observable 
for only a short time. The rarer binary neutron stars, like PSR 1913-1-16, do 
have a long-lived recycled pulsar, which more than offsets their lower formation 
rate and makes them dominate the observed population. 

We do observe high-mass black holes in CygX-1 and in soft X-ray transients. 
In the former, the black hole is of > 10 Mq [4]. The companion 0-star is near 
its Roche Lobe, and its wind is continuously feeding the black hole, which 
shines through X-ray emission. In addition to CygX-1, high-mass black holes 
are seen in the LMC in LMCX-3 and perhaps LMCX-1. Much more copious 
are the transient sources, with black holes of mass Mbh ~ 7 M©, most of 
which flare up only occasionally with long quiescent times between flare ups. 
Wijers [5] estimated ~ 3000 of these in the Galaxy. That is, these are the 
numbers that are presently operative. Remarkable about the transient sources 
with unevolved donors is that the main sequence star is K- or M-star, less 
massive than our sun. Brown, Lee, & Bethe [6] explain this in terms such 
that higher-mass donors can also participate in the formation of such binaries 
containing a high-mass black hole, but will end up in the evolution further 
away from the black hole so that they can pour matter on the latter only 
when they evolve in subgiant or giant stage. Thus, there are a large factor 
estimated to be ~ 22 more of those binaries which will not be seen until the 
main sequence star evolves [6] . The mechanism describing the evolution of the 
transient sources required the massive progenitor of the black hole to carry 
out core helium burning as if it were a single star; i.e., before having its H 
envelope removed in RLOF by its main sequence companion. An interval of 
~ 20 — 35 Mq ZAMS was estimated for the progenitors of the high-mass black 
hole. Consequently, this same interval of single stars, not in binary, would be 
expected to end up as high-mass black holes. In the formation of these high- 
mass black holes, most of the helium envelope of the progenitor must drop into 
the black hole in order to form their high mass, so little matter is returned to 
the Galaxy. 

This brings us to the intriguing matter of SN 1987A which we believe did go 
into a black hole, but after the supernova explosion which did return matter 
to the Galaxy. The progenitor of SN 1987A was known to have ZAMS mass 
~ 18 Mq. This leads us to the interesting concept of low-mass black holes with 
delayed explosion, which result from the ZAMS mass range ~ 18 — 20 Mq, 
although the precise interval is uncertain. The delayed explosion mechanism 
has been elucidated by Prakash et al. [7]. The range of ZAMS masses of single 
stars in which neutron stars are formed is thus only ~ 10 — 18 Mq. 

The absence of matter being returned to the Galaxy in the ZAMS mass range 
~ 20 — 35 Mq impacts on nucleosynthesis, especially in the amount of oxygen 
produced. Bethe & Brown [8] suggested that matter was again returned to 



the Galaxy by stars in the ZAMS range ~ 35 — 80 Mq. In this case, the 
progenitor was stripped of H envelope in an LBV phase, and the naked He 
star was suggested to evolve to a low-mass black hole, with return of matter 
to the galaxy before its formation in a delayed explosion, or to a neutron star. 
Thus, elements like oxygen were produced in a bimodal distribution of ZAMS 
masses M < 20 M© and 35 M© < M < 80 Mq. 

The Bethe & Brown [8] suggestion was based on naked He stars evolved by 
Woosley, Langer, & Weaver [9] who used a too-large wind loss rate for He stars. 
Wellstein & Langer [10] have evolved naked He stars with lower rates, in which 
case the final He envelope is somewhat larger. However, the central carbon 
abundance following core He burning is high ~ 33%. With this abundance, 
the stars will not skip the convective carbon burning stage in their evolution, 
and according to the arguments of Brown, Lee, & Bethe [6] would still be 
expected to end up as low-mass compact objects, in which case matter would 
be returned to the Galaxy. This matter will not, however, be settled until the 
CO cores evolved with lowered He-star wind loss rates by Wellstein & Langer 
have been burned further up to the Fe core stage, so the Bethe & Brown [8] 
bimodal mass region for nucleosynthesis should be viewed as provisional. 

In sect. 2, we discuss the maximum mass of neutron stars and the processes 
that determine which range of initial stellar masses gives rise to what compact 
object, and how mass loss in naked helium stars changes those ranges. Then 
we describe the Bethe & Brown [1] scenario for the evolution of massive bi- 
nary stars, and especially their treatment of common-envelope evolution and 
hypercritical accretion (sect. 3). We then discuss a few specific objects sepa- 
rately, first binary neutron stars (sect. 4), then CygX-1 and its ilk (sect. 5) and 
the black-hole transients (sect. 6). Then we comment briefly on how our re- 
sults would affect the evolution of low-mass X-ray binaries with neutron stars 
(sect. 7) and summarize our conclusions (sect. 8). The discussion of CygX-3 
and the possible implications of neutron-star, black-hole binaries for gravity 
waves and gamma-ray bursts are in appendices 1-3. 



2 The Compact Star 



Thorsson, Prakash, & Lattimer [11] and Brown & Bethe [12] have studied the 
compact core after the collapse of a supernova, assuming reasonable interac- 
tions between hadrons. Initially, the core consists of neutrons, protons and 
electrons and a few neutrinos. It has been called a proto-neutron star. It is 
stabilized against gravity by the pressure of the Fermi gases of nucleons and 
leptons, provided its mass is less than a limiting mass Mpc (proto-compact) 
of ~ 1.8M0. 



If the assembled core mass is greater than Mpc there is no stabihty and no 
bounce; the core collapses immediately into a black hole. It is reasonable to 
take the core mass to be equal to the mass of the Fe core in the pre-supernova, 
and we shall make this assumption, although small corrections for fallback in 
the later supernova explosion can be made as in Brown, Weingartner & Wijers 
[13]. If the center collapses into a black hole, the outer part of the star has no 
support (other than centrifugal force from angular momentum) and will also 
collapse. 

If the mass of the core is less than Mpc, the electrons will be captured by 
protons 



p + e ^ n + V (1) 

and the neutrinos will diffuse out of the core. This process takes of order of 
10 seconds, as has been shown by the duration of the neutrino signal from SN 
1987A. The result is a neutron star, with a small concentration of protons and 
electrons. The Fermi pressures of the core are chiefly from the nucleons, with 
small correction from the electrons. On the other hand the nucleon energy is 
increased by the symmetry energy; i.e., by the fact that we now have nearly 
pure neutrons instead of an approximately equal number of neutrons and 
protons. Thorsson et al. [11] have calculated that the maximum mass of the 
neutron star Mns is still about l.SM©; i.e., the symmetry energy compensates 
the loss of the Fermi energy of the leptons. Corrections for thermal pressure 
are small [14]. 

The important fact is that the ten seconds of neutrino diffusion from the core 
give ample time for the development of a shock which expels most of the mass 
of the progenitor star. 

But this is not the end of the story. The neutrons can convert into protons 
plus K^ mesons. 



n^p + R-. (2) 

This is short-hand for the more complicated interaction iV + e" ^ N' + K~ + v 
where iV is a nucleon. The neutrinos leave the star. The times are sufficiently 
long that chemical equilibrium is assured. Since the density at the center of 
the neutron star is very high, the energy of the K~ is very low, as confirmed 
by Li, Lee, & Brown [15] using experimental data. By this conversion the 
nucleons can again become half neutrons and half protons, thereby saving the 
symmetry energy needed for pure neutron matter. The K~ , which are bosons, 
will condense, saving the kinetic energy of the electrons they replace. The 
reaction eq. (2) will be slow, since it is preceded by 



e- ^ R- + u (3) 

(with the reaction eq. (2) following) as it becomes energetically advantageous 
to replace the fermionic electrons by the bosonic K^^s at higher densities. 
Initially the neutrino states in the neutron star are filled up to the neutrino 
chemical potential with trapped neutrinos, and it takes some seconds for them 
to leave the star. These must leave before new neutrinos can be formed from 
the process eq. (3). Thorsson et al. [11] have calculated that the maximum 
mass of a star in which reaction eq. (2) has gone to completion is 

Mnp ~ 1.5 Mq , (4) 

where the lower suffix NP denotes their nearly equal content of neutrons and 
protons, although we continue to use the usual name "neutron star" . This is 
the maximum mass of neutron stars, which is to be compared with the masses 
determined in binaries. The masses of 19 neutron star masses determined in 
radio pulsars [16] are consistent with this maximum mass. 

The core mass Mc formed by the collapse of a supernova must therefore be 
compared to the two hmiting masses, Mpc and Mnp. If 

(I) Mc > Mpc (5) 
we get a high mass black hole. If 

(II) Mpc > Mc > Mnp (6) 
we get a low- mass black hole, of mass Mc. Only if 

(III) Mc < Mnp (7) 

do we get a neutron (more precisely, "nucleon") star from the SN. Only in this 
case can we observe a pulsar. In cases (II) and (III) we can see a supernova 
display. In case (I) we receive only initial neutrinos from electrons captured in 
the collapse before Mc becomes greater than Mpc but no light would reach 
us. (Except perhaps if the new black hole rotates rapidly enough to power an 
explosion, a mechanism proposed by MacFadyen and Woosley [93] for gamma- 
ray bursts.) 

Woosley, Langer, & Weaver [17] evolve massive stars with mass loss. For stars 
in the ZAMS mass range ~ 20 — 30 Mq, mass loss is relatively unimportant 
and since Mpc ^1.8 Mq for this range, we find from the earlier calculation of 
Woosley & Weaver [18] that most of the single stars in this range will go into 
high-mass black holes. Evolution of these stars in binaries is another matter. 



Timmes, Woosley, & Weaver [19], Brown, Weingartner, & Wijers [13], and 
Wellstein & Langer [10] find that substantially smaller core masses result if the 
hydrogen envelope is taken off in RLOF so that the helium star is naked when 
it burns. Thus, stars of ZAMS masses ~ 20 — 35 Mq in such binaries evolve 
into low-mass compact cores, black hole or neutron star. Woosley, Langer, 
& Weaver [17] used helium-star wind loss rates which were too high by a 
factor ~ 2 — 3, but lower wind losses give only slightly larger He cores in the 
ZAMS mass range ~ 20 — 35 Mq [10] so our above conclusion is unlikely to 
be reversed. 

On the other hand, the fate of single stars in the ZAMS mass range ~ 35 — 
80 Mq is uncertain. In the published Woosley, Langer, & Weaver [17] work 
with too high mass loss rate, so much matter is blown away, first in LBV stage 
and later in W.-R. stage that low-mass compact objects, black-hole or neutron- 
star, result [13]. Bethe & Brown [8] attribute this to the fact that convective 
carbon burning is not skipped in these stars. In this stage a lot of entropy can 
be removed by vv emission, so that a low-entropy, and therefore small, core 
results. In this range, Wellstein & Langer [10] find central ^^C abundances of 
33 — 35% following He core burning, more than double the ~ 15% required for 
convective carbon core burning. Therefore, we believe that this range of stars 
will still go into low-mass compact objects, even though their final He cores 
are substantially larger because of the lower, more correct. He-star wind mass 
loss rates used by Wellstein & Langer [10]. However, this problem cannot be 
considered as settled until the Wellstein & Langer CO cores are burned up 
to the Fe core stage. We will therefore not discuss the evolution of CygX-1 
like objects, high-mass black holes accompanied by sufficiently massive giant 
companion so that they shine continuously in X-rays. It is not clear to us 
whether LMCX-3, with a high-mass black hole and a B-star companion of 
roughly equal mass, has a history more like CygX-1 or like the transient 
black-hole binaries which we discuss below. 

Bethe & Brown [8] took 80 Mq as the lower mass limit for high-mass black hole 
formation in binaries which experience RLOF; i.e., in those for which helium 
core burning proceeds in a naked helium star. Because of our above discussion, 
we believe this mass limit may be too high, so that the contributions from 
high-mass black-hole, neutron-star binaries were, if anything, underestimated 
in their work. However, we will not know until the CO cores obtained with 
better He-star mass loss rates are evolved further. 



3 Evolution of Binary Compact Objects 



We summarize the Bethe & Brown [1] evolution of binary compact objects, 
paying special attention to their common envelope evolution. In particular, we 
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shall show that their schematic evolution should be applicable to donors with 
deeply convective envelopes, whereas for non-convective or semiconvective en- 
velopes, such as encountered in the evolution of low-mass X-ray binaries, their 
common envelope evolution would not be expected to apply. 

We call the star that is initially heavier star A, the other star B. We denote 
initial masses by subscript i, so we have masses Ma,j, Mb,^. We denote their 
ratio by q; thus 

g = Mb,,/Ma,. < 1 . (8) 

Following Portegies Zwart & Yungelson [20], we assume that q is distributed 
uniformly between and 1. Likewise, we also assume that In a is uniformly 
distributed, where a is the semi-major axis of their orbit. 

However, we assume different limits for a than Portegies Zwart & Yungelson 
[20] . Initially both stars are massive main sequence stars, with radius at least 
3 Rq, so a > 6 Rq = 4 X 10^ km. At the other end of the scale, we require 
a < 4 X 10^ km. We assume that 50% of all stars are binaries with separations 
in this range (stars in wider binaries would evolve as if they were single). Then 
the fraction of binaries in a given interval of In a is 

d0 = ci(lna)/7. (9) 

We assume that a star needs an initial mass of 

M > M, = 10 Mq (10) 

in order to go supernova. Therefore, if a is the total rate of SNe, the rate of 
SNe in mass interval dM is given by 



( M \ "" dM 

da = an --— — — - (11) 

\10 MqJ M ^ ' 

where we have used a power-law initial mass function with n = 1.5 (close to 
the Salpeter value n = 1.35). The birth rate of supernova systems was taken 
to be 



a = 0.02yr-^ (12) 

in the Galaxy. By a supernova system we mean a single star that goes super- 
nova (i.e., has Mzams > 10 Mq) or a close binary containing at least one such 
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star (close here means within the separation range mentioned above). Bethe 
& Brown [1] find that if the primary is massive enough to go supernova, then 
there is an ~ 50% chance for the secondary to also go supernova. This was 
calculated for a distributions fiat in g = M-Q^i/Mp^i. Therefore, the supernova 
rate in our notation would be 1.25a = 0.025 yr^^. 

Using the Cordes & Chernoff [21] distribution of kick velocities, 43% of the 
binaries were found to survive the first explosion. Thus, at this stage, we are 
left with a birth rate of 



R = 0.02 X - X - X 0.43 ^ 2 x 10^^ per yr (13) 

for the formation of binaries consisting of a neutron star with a companion 
massive enough to go supernova (M > 10 Mq). The life time of such systems 
is the companion life time of ~ 10^ yr, but star A will be a pulsar for only 
~ 5 x 10® yr because it will spin down electromagnetically until it is no longer 
observable. From these numbers we estimate the number of such systems to 
be ~ 10"^ in the Galaxy. 

Since the pulsar is unrecycled, the expected number should be compared with 
the detected population of active radio pulsars in the galaxy, about 10^. This 
number should be multiplied by a factor of 1/2 for binarity, a further factor of 
1/2 for a binary in which both stars can go supernova and the 0.43 for survival 
of the first explosion. This would leave the large number ~ 10^ if pulsars with 
massive companions were as easily detected as single pulsars. In fact, only 2 
are observed; PSR1259— 63 with a Be-star companion and PSR0045— 73 with 
a B-star companion. Stellar winds interfere with the radio pulses from these 
binaries, obscuring the narrower ones. Doppler shifts also make these difficult 
to observe. Nevertheless, the factor necessary to reduce their observability is 
large. We return to the subject later. 

At this stage we have an ~ 1.4 Mq neutron star with O or B-star companion. 
We take the latter to have mass ~ 15 Mq. The giant has a He core containing 
some 30% of its mass, surrounded by an envelope consisting mainly of H. 
We take the envelope to be deeply convective,Q so the entropy is constant. 
The particles, nuclei and electrons, are nonrelativistic and thus have 7 = 5/3. 
Therefore, the envelope forms a polytrope of index n = 3/2. Applegate [26] 
shows that the binding energy of the envelope is 



^ The assumption that the envelope is deeply convective is essential for our later 
treatment of common envelope evolution with hypercritical accretion. Recent devel- 
opments with non-convective or semiconvective donors show that the accretion rate 
is also highly super-Eddington, but still significantly less [22-24]. For very massive 
donors the rate is always highly super-Eddington. 
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E ^ O.GGM^R-' (14) 

where R is the outer radius. In this formula the binding energy is decreased 
50% by the kinetic energy, E containing both effects. 

The major difference of the Bethe & Brown calculations and of case H of 
Portegies Zwart & Yungelson [20] compared with other work is the use of hy- 
percritical accretion. In a series of papers Chevalier [27,28] showed that once M 
exceeded ~ 10^ M-^^^, the photons were carried inwards in the adiabatic inflow, 
onto the neutron star. The surface of the latter was heated sufficiently that 
energy could be carried off by neutrino pairs. Brown [3] reproduced Cheva- 
lier's results in analytical form. The idea has a much longer history: Colgate 
[73] showed already in 1971 that if neutrinos carry off the bulk of the energy, 
accretion can proceed at a much greater rate than Eddington. In 1972 Zel- 
dovich et al. [74], before the introduction of common envelope evolution, used 
hypercritical accretion of a cloud onto a neutron star. Bisnovatyi-Kogan & 
Lamzin [75] and Chevalier [27] pointed out that during the common envelope 
phase of binary evolution, photons would be trapped and accretion could oc- 
cur at much higher rates, and that neutron stars that go through this phase 
generally will go into black holes. 

We begin by considering the work done by the neutron star on the envelope 
matter that it accretes. This will turn out to be only a fraction of the total 
work, the rest coming from the production of the wake, but it illustrates simply 
our procedure. Taking the neutron star to be at rest, the envelope matter is 
incident on it with the Keplerian velocity v. The rate of accretion is given by 
the Bondi-Hoyle-Lyttleton theory 



^ = -F.flL (15) 

where p is the density of the B material, v is its velocity relative to the neutron 
star A, and R^c is the accretion radius 



i?ac = 2GMav"^. (16) 

The rate of change of momentum P is 



dP dMA ,^^. 

11 = ""^r^ ('^) 

the matter being brought to rest on the neutron star, and this is equal to the 
force F. Consequently, the rate at which the neutron star does work in the 
material is 
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h = hv = V 



dt 



Inclusion of the work done in creating the wake involves numerical calculations 
[29-31] with the result that the coefficient of the right-hand side of eq. (18) is 
changed; i.e, 



Cd\ 2dMA 



with Q ~ 6 — 8 for our supersonic flow. It is, in fact, very important that the 
wake plays such a large role, in that it's the fact that q/2 > 1 (We consider 
q/2 to be 3> 1) that makes our later common envelope evolution strongly 
nonconservative, the proportion of the total H-envelope mass accreted onto 
the neutron star being relatively small. 

In eq. (19) v^ is the velocity of the B (giant) material relative to A, the neutron 
star. This is given by 



v^ = G{Ma + Mb) a-\ (20) 

The interaction energy of A and B is 

E = -GMAMBa-\ (21) 

Since we know Mb,^ and Mb,/, the initial mass of B and the mass of its He 
core, our unknown is a/. We can obtain it by considering 

Y = M^a-^ (22) 

as one variable. Ma as the other. Differentiating eq. (21) we have 



E = ]f (MpJ + Mpj) (23) 

whereas combining eqs. (19) and (20) and neglecting Ma with respect to Mb, 
we have 



E = g{^^^Y Ma. (24) 

Thus, eqs. (23) and (24) are equal, so we have 
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Ma _ 1 F 
which can be integrated to give 



Ma oc ri/(^''-i) = yi/s (26) 

where we have chosen c^ = 6 [29]. The final energy is then 



Ef = -GM^,YJ-l) . (27) 



2 '\YJ 



The binding energy Ef of star A to star B serves to expel the envelope of 
star B, whose initial binding energy is given by eq. (14). Mass transfer begins 
at the Roche Lobe which lies at ~ 0.6 a^ for the masses involved. However, 
star B expands rapidly in red giant stage before the mass transfer can be 
completed. To keep the numbers easy to compare with Bethe & Brown [1], 
we use their approximation of starting spiral-in when the giant's radius equals 
the orbital separation rather than the Roche-lobe radius. Since for the large 
mass ratios considered here, Ri^/a ~ 0.5 for the giant, this implies we require 
Ef of eq. (27) to be about twice the binding energy (eq. 14), i.e. 

E, = Mg<1 = 1,2 G^. (28) 

(We set the common-envelope efficiency, a, to 0.5.) The ejected material of B 
is, therefore, released with roughly the thermal energy it had in the envelope; 
in other words, the thermal energy content of the star is not used to help expel 
it. Inserting eq. (28) into eq. (27) yields 



Star A is initially a neutron star, Ma,j = 1.4 Mq. For star B we assume 
MB,i = 15 Mq. Then eq. (29) yields 



^ = 15. (30) 

We use this to find the result of accretion, with the help of eq. (26), 
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^ = 1-73 (31) 

MA,i ^ ^ 



or 



Maj = 2.4 Mq. (32) 

This is well above any of the modern limits for neutron star masses, so we find 
that the neutron star has gone into a black hole. 

Our conclusion is, then, that in the standard scenario for evolving binary neu- 
tron stars, if the giant is deeply convective, accretion in the common envelope 
phase will convert the neutron star into a black hole. 

Star B, by losing its envelope, becomes a He star. We estimate that 

0.3. (33) 



Mb,/ 



Mb,. 
The size of the orbit is determined by eq. (22) 



^ = ^^^ = 50. (34) 

The final distance between the stars a/ should not be less than about 10^^ cm, 
so that the He star (mass Mb,/) fits within its Roche lobe next to the black 
hole of mass Ma,/. Bethe & Brown [1] showed that if the black hole and the 
neutron star resulting from the explosion of star B are to merge in a Hubble 
time, then a/ < 3.8 x 10^^ (for circular orbits; correction for eccentricity will 
be given later). Therefore the initial distance of the two stars, after the first 
mass exchange and the first supernova should be 

0.5 X lO^''^ cm < Oi < 1.9 x 10^^ cm (35) 



If the initial distribution of distances is da/la, the probability of finding a 
between the limits of eq. (35) is 

P = 18%. (36) 

As noted earlier, 43% of the binaries survive the first explosion, so the com- 
bined probability is now 

P = 8% (37) 



16 



for the survivors falling in the logarithmic interval in which they survive coa- 
lescence, but are narrow enough to merge in a Hubble time. Our final result, 
following from a birth rate of 10~^ binaries per year in which one star goes 
supernova, half of which have both stars going supernova, is 



R = IQ-^ X 0.5 X 0.08 X 0.5 = 2 X 10"^ yr~' (38) 

in the Galaxy. The final factor of 0.5 is the survival rate of the He-star, neutron 
star binary, calculated by Monte Carlo methods. Bethe & Brown [1] quoted 
10"*^ yr"^, or half of this rate, in order to take into account some effects not 
considered by them in which the binary disappeared (e.g., Portegies Zwart 
and Verbunt [57]). 

Our final rate is, then. 



R = 10- V"^ galaxy-^ (39) 

Using our supernova rate of 0.025 per year, which includes the case where 
both stars in the binary go supernova, we can convert this birth rate to 0.004 
per supernova for comparison with other work. Portegies Zwart & Yungelson 
[20] in their case H, which included hypercritical accretion, got 0.0036 per 
supernova, within 10% of our value. Thus, the chief difference between our 
result in eq. (39) and the R = 5.3 x 10^^ of these authors is due to the 
different assumed SN rate. 

In our above estimates we have assumed the second neutron star to be formed 
to have a circular orbit of the same a as its He-star progenitor. However, 
eccentricity in its orbit leads to a value of aj substantially larger than the 3.8 x 
10^^ cm used above as the maximum separation for merger. In general, most of 
the final binaries will have e > 0.5, with a heavy peak in the distribution close 
to e = 1. The rise occurs because preservation of the binary in the explosion 
is substantially greater if the kick velocity is opposite to the orbital velocity 
before explosion. In this case the eccentricity e is large. The most favorable 
situation is when the orbital and kick velocities are equal in magnitude. (See 
the figures in Wettig and Brown [32].) Eggleton [33] has kindly furnished us 
with a useful interpolation formula for the increase. The factor by which to 
multiply the time for merger in circular orbits, is 



Z(e) f« (1 - e')3-^^9-°-243^-°-°^^'= . (40) 

This formula is accurate to about 1% for e < 0.99. Thus, if the initial eccen- 
tricity is 0.7, the time to shrink the orbit to zero is about 10% of the time 
required if the initial eccentricity were zero for the same initial period. The 
maximum aj = 3.8 x 10^^ cm for circular orbits would be increased by the 
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fourth root of the decrease in time; i.e., up to 6.8 x 10^^ cm for this eccentric- 
ity. The maximum Oj in eq. (35) would go up to 3.4 x 10^^ cm, increasing the 
favorable logarithmic interval by ~ 40%. We have not introduced this correc- 
tion because it is of the same general size as the uncertainty in the supernova 
rate. However, this correction gives us some comfort that our final numbers 
are not unreasonably large. 

If we produce an order of magnitude more low-mass black hole, neutron star 
binaries than binary neutron stars, the obvious question is why we have not 
seen any. The neutron star in this object is "fresh" (unrecycled) so it would 
spin down into the graveyard of neutron stars in ~ 5 x 10^ yr. The two 
relativistic binary pulsars we do see 1913+16 and 1534+12 have been recycled, 
have magnetic fields B ~ 10^'^ G, two orders of magnitude less than a fresh 
pulsar, and will therefore be seen for about 100 times longer than an unrecycled 
neutron star. So even with a ten times higher birth rate, we should see ten 
times fewer LBH-NS binaries than NS-NS binaries. Furthermore, the binary 
with black hole will have a somewhat higher mass, therefore greater Doppler 
shift, and therefore be harder to detect. In view of the above, it is reasonable 
that our low-mass black-hole, neutron-star binaries have not been observed, 
but they should be actively looked for. 

We should also calculate the rate of coalescences of the black hole with the 
He star. These have been suggested by Fryer & Woosley [34] as candidate 
progenitors for the long time gamma ray bursters. Note that they will occur 
for a range of 0.04 x 10^^ cm < a^ < 0.5 x 10^^ cm, a logarithmic interval double 
that of eq. (35). Thus, the black-hole. He-star coalescence has a probability 

P = 36%. (41) 

Furthermore, this situation does not have the 50% disruption in the final 
explosion, so the black-hole. He-star coalescences occur with a total rate of 4 
times that of the black-hole, neutron-star mergers. 

There has been much discussion in the literature of the difficulties in com- 
mon envelope evolution. We believe our model of deeply convective giants and 
hypercritical accretion offers an ideal case. Of course, the initiation of the 
common envelope evolution requires some attention, but it can be modeled in 
a realistic way [35]. As the giant evolves across its Roche lobe, the compact 
object creates a tidal bulge in the giant envelope, which follows the com- 
pact object, torquing it in. As the convective giant loses mass, the envelope 
expands in order to keep entropy constant. In Bondi-Hoyle-Lyttleton accre- 
tion, a density poo ~ 10^^^ g cm^'^ is sufficient with wind velocities ~ 1000 
km s~^ in order to give accretion at the Eddington rate. Thus to achieve 
M ~ lO^MEdd ~ 1 MqJt"-^ we need p ~ 10"^ g cm~^ which is found at 0.9 R, 
where R is the radius of the giant. At this rate of accretion, angular momen- 



turn, etc, are hardly able to impede it appreciably. The total mass accreted 
onto the compact object is ~ 1 Mq, so the common envelope evolution has 
dynamical time of years. As noted earlier, it is non-conservative. 



4 Evolution of Binary Neutron Stars 



Since the standard scenario of evolution of binary compact objects ends up 
with low-mass black-hole, neutron-star binaries, another way must be found 
to evolve neutron star binaries. In the double He-star scenario was suggested 
by Brown [3] and developed further by Wettig & Brown [32] the neutron star 
avoids going through common envelope with a companion star. In this way 
the neutron star can avoid being converted into a black hole by accretion. For 
two giants to burn He at the same time, they must be within ~ 5% of each 
other in mass, the helium burning time being ~ 10% of the main sequence 
life time, and stellar evolution time going roughly with the inverse square 
of the mass. With a flat mass ratio distribution, this happens in 5% of all 
cases, making the ratio of NS-NS to NS-LBH binaries 1:20. However, when 
the primary becomes an LBH, only half the secondaries will be massive enough 
to form a NS, whereas for the very close mass values of the double-He scenario 
this factor 2 loss does not occur. Thus, binary neutron stars should be formed 
10% as often as low-mass black-hole, neutron-star binaries. This 10% is nearly 
model independent because everything else roughly scales. 

The scenario goes as in Wettig & Brown [32]. The primary 0-star evolves 
transferring its H-envelope to the companion. Often, this would lead to 'reju- 
venation' of the secondary, i.e. its evolution would restart also from the ZAMS 
with the now higher total mass, and it would make a much heavier core. How- 
ever, here the core of the secondary has evolved almost as far as the primary's 
core, so the core molecular weight is much higher then that of the envelope. 
This prevents convection in the core from extending into the new envelope to 
make the bigger core, so no rejuvenation takes place [36]. Since g' ~ 1, the first 
mass transfer is nearly conservative. The second is not, so the two He-cores 
then share a common H envelope, which they expel, while dropping to a lower 
final separation aj. 

Following the explosion of the first He star, the companion He-star pours wind 
matter onto the pulsar, bringing the magnetic field down and spinning it up 
[3,32]. The end result is two neutron stars of very nearly equal mass, although 
wind accretion can change the mass two or three percent. 

The above scenario ends for He-star masses greater than 4 or 5 Mq, corre- 
sponding to ZAMS masses greater than ~ 16 or 18 Mq. However, less massive 
He stars evolve in the He shell-burning stage, and a further mass transfer (Case 
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C) can take place. The transfer of He to the pulsar can again bring about a 
black hole, which Brown [37] very roughly estimates to occur in ~ 50% of the 
double neutron star binaries. This is roughly consistent with results of Fryer & 
Kalogera [38]. Taking a rate oi R = 10~"^ per year per galaxy for the low-mass 
black-hole, neutron-star binaries, we thus arrive at a birth rate of 



i? ~ 5 X 10 ^ per year per galaxy (42) 

for binary neutron-star formation. However, the black holes formed in the He 
shell burning evolution will not have accreted much mass and will have about 
the same chirp mass as binary neutron stars (see below) for gravitational 
merging. 

Our best guess values eqs. (39) and (42) thus give an ~ 20 to 1 ratio for forma- 
tion of low-mass black-hole, neutron-star binaries to binary neutron stars. The 
former are better progenitors for gravitational waves from mergers because of 
their higher masses and they have many advantages as progenitors of gamma 
ray bursters [39]. Note that our estimated rate of i? = 5 x 10^^ per galaxy per 
year for binary neutron star formation is consistent with the empirical rates 
discussed in our introduction. 



5 High-Mass Black-Hole O/B-star Binaries 



We will be brief in our review of these, because we believe the evolution of these 
objects such as CygX-1, LMCX-1 and LMCX-3 to be less well understood 
than the low-mass black-hole, neutron-star binaries. Evolutionary calculations 
now proceeding by Alexander Heger, using the CO cores evolved by Wellstein 
& Langer [10] should clarify this situation substantially. 

Bethe & Brown [8] arrived at a limit of ZAMS mass 80 Mq for stars in binaries 
to go into high-mass black holes (unless Case C mass transfer takes place 
as we discuss in our next section). This limiting mass is much higher than 
other workers have used. It was based on calculations of Woosley, Langer, & 
Weaver [17] and was so high because of very high mass loss rates used by these 
authors. With more correct lower rates the limiting mass may come down, so 
the Bethe & Brown evolution should be viewed as giving a lower limits to 
the number of high-mass black-hole, O/B-star binaries. Their estimated birth 
rate of about 3 x 10"^ per galaxy per year does agree reasonably well with the 
fact that only one such system is known in the Galaxy. However, since even 
with a twice larger separation the accretion rate of the black hole from the 
fast wind of the O star becomes small, it is possible that substantially more 
systems with somewhat wider orbits exist undetected, and that CygX-1 is the 
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only one presently in the (very short) phase of incipient Roche lobe overflow 
when it is bright. Bethe & Brown [8] found this narrowness of the CygX- 
1 orbit (17 Rq according to Herrero et al. [4]) to be puzzling: the massive 
stars in the progenitor binary initially had to flt within their Roche Lobes, 
therefore a separation of at least double the current 17 Rq was needed. And 
most evolutionary effects from then on, such as wind mass loss or supernova- 
like mass loss, would tend to widen the orbit. Of course, the orbit could be 
narrowed in Case A mass transfer (i.e. during the main sequence) since the 
progenitor of the black hole was more massive than the present donor, but it 
could not become so narrow that the present donor fllled its Roche lobe, and 
would widen again once the mass ratio became reversed and widen further 
due to wind loss after the whole primary envelope was lost. 

In any case, a binary as narrow as CygX-1 would coalesce in the common 
envelope evolution once the 0-star companion of the massive black-hole goes 
into red giant phase, according to the Bethe & Brown [8] estimates. Since the 
black hole in CygX-1 has mass > 10 Mq and is probably the most massive 
black hole in a binary observed in the Galaxy, in the Fryer & Woosley [34] 
model where the black hole "eats" the W.-R. companion, such a coalescence 
should produce the most energetic long-lasting gamma ray burster. We are 
unable to evaluate the probability of CygX-1 like objects merging following 
common envelope evolution because we have been unable to understand why 
CygX-1, before common envelope evolution, is so narrow. The LBV, RSG, 
and WNL stages of W.-R. development are not quantitatively understood. 

After the main sequence star in a Cyg X-l-like object explodes and becomes 
a neutron star, according to Bethe & Brown [8] the binary will eventually 
merge. They estimated the contribution to the merger rate of these systems 
to be (4 — 6) X 10^® yr^^ galaxy^-*^, however with considerable uncertainty due 
to the fact that the evolution of CygX-1 itself is uncertain. Lowering the mass 
limit for black-hole formation by having lower mass loss rates would increase 
this number (e.g. a limit of 40 Mq would increase the merger rate by a factor 
5). 



6 The Formation of High-Mass Black Holes in Low-Mass X-ray 
Binaries 



6. 1 General 



Crucial to our discussion here is the fact that single stars evolve very differently 
from stars in binaries that lose their H-envelope [6,10,13,19] either on the main 
sequence (Case A) or in the giant phase (Case B). However, stars that transfer 
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mass or lose mass after core He burning (Case C) evolve, for our purposes, 
as single stars, because the He core is then exposed too close to its death for 
wind mass loss to significantly alter its fate. Single stars above a ZAMS mass 
of about 20 Mq skip convective carbon burning following core He burning, 
with the result, as we shall explain, that their Fe cores are substantially more 
massive than stars in binaries, in which H-envelope has been transferred or 
lifted off before He core burning. These latter "naked" He stars burn ^^C 
convectively, and end up with relatively small Fe cores. The reason that they 
do this has to do chiefly with the large mass loss rates of the "naked" He 
cores, which behave like W.-R.'s. Unfortunately, in calculation until recently, 
substantially too large mass loss rates were used, so we cannot pin limits 
down quantitatively. In this section we will deal with the ZAMS mass range 
~ 20 — 35 Mq, in which it is clear that many, if not most, of the single stars 
go into high-mass black holes, whereas stars in binaries which burn "naked" 
He cores go into low-mass compact objects. In this region of ZAMS masses 
the use of too-high He-star mass loss rates does not cause large effects [6]. 

The convective carbon burning phase (when it occurs) is extremely important 
in pre-supernova evolution, because this is the first phase in which a large 
amount of entropy can be carried off in z/z/-pair emission, especially if this 
phase is of long duration. The reaction in which carbon burns is ^^C(a, 7)^^0 
(other reactions like C + C would require excessive temperatures). The cross 
section of ^^C(q;,7)^^0 is still not accurately determined; the lower this cross 
section the higher the temperature of the ^^C burning, and therefore the more 
intense the vv emission. With the relatively low ^^C(a, 7)^^0 rates determined 
both directly from nuclear reactions and from nucleosynthesis by Weaver & 
Woosley [42], the entropy carried off during ^^C burning in the stars of ZAMS 
mass ~ 10 — 20 Mq is substantial. The result is rather low-mass Fe cores for 
these stars, which can evolve into neutron stars. Note that in the literature 
earlier than Weaver & Woosley [42] often large ^^C(q;,7)^^0 rates were used, 
so that the ^^C was converted into oxygen and the convective burning did not 
have time to be effective. Thus its role was not widely appreciated. 

Of particular importance is the ZAMS mass at which the convective carbon 
burning is skipped. In the Woosley & Weaver [18] calculations this occurs at 
ZAMS mass IQM© but with a slightly lower ^^C(a, 7)^^0 rate it might come at 
2OM0 or higher [37]. As the progenitor mass increases, it follows from general 
polytropic arguments that the entropy at a given burning stage increases. At 
the higher entropies of the more massive stars the density at which burning 
occurs is lower, because the temperature is almost fixed for a given fuel. Lower 
densities decrease the rate of the triple-a process which produces ^^C relative 
to the two-body ^^C(a, 7)^^0 which produces oxygen. Therefore, at the higher 
entropies in the more massive stars the ratio of ^^C to ^^O at the end of He 
burning is lower. The star skips the long convective carbon burning and goes 
on to the much shorter oxygen burning. Oxygen burning goes via ^^O +^'^ O 



22 



2.4 



2.2 - 



2.0 - 



^ 1-8 

CO 



o 1.6 



1.4 



1.2 
10.0 



/ 


- 






- 


■J 


PSRB1913+16 


- 



15.0 20.0 

ZAMS Mass I 



25.0 



30.0 



Fig. 1. Compact core masses resulting from the evolution of single stars, Case B 
of solar metallicity of Woosley & Weaver [18]. The horizontal dashed lines indicate 
the mass of the heaviest known well-measured pulsar [16], the maximum mass of 
a neutron star, and our estimate of Mpc (proto-compact), the maximum compact 
core mass for which matter can be returned to the galaxy. 

giving various products, at very much higher temperature than C(q;,7) and 
much faster. Since neutrino cooling during the long carbon-burning phase gets 
rid of a lot of entropy of the core, skipping this phase leaves the core entropy 
higher and the final Chandrasekhar core fatter. 

In Fig. 1 the large jump in compact object mass in single stars at ZAMS mass 
~ 19 Mq is clearly seen. From our discussion in Section 2 we see that this is 
just at the point where the Fe core mass goes above the proto-compact mass 
of ~ I.SMq and, therefore, above this mass one would expect single stars to go 
into high-mass black holes. Arguments have been given that SN 1987A with 
progenitor ZAMS mass of ~ 18 Mq evolved into a low-mass black hole [12]. 
We believe from our above arguments and Fig. 1 that just above the ZAMS 
mass of ~ 20 Mq, single stars go into high- mass black holes without return 
of matter to the Galaxy. Thus, the region of masses for low-mass black hole 
formation in single stars is narrow, say ~ 18 — 20 Mq (although we believe it 
to be much larger in binaries). 

Thus far our discussion has been chiefly about single stars, in which the He 
burns "clothed" by a hydrogen envelope. In this case the convective helium 
core grows in stars as time passes. In the "naked" He cores, in which the H 
envelope has been lifted off in RLOF or driven off by wind either before or 
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early in the He burning the temperature and the entropy will be slightly lower, 
because the insulating layer is gone, so it is not surprising that their carbon 
abundance is large. Furthermore, the core mass continually decreases because 
of mass loss by wind. In fact, even for the naked 20 Mq He core, corresponding 
to ZAMS mass 45 M©, the central carbon abundance was ~ 33% at the end 
of He core burning [18] whereas only ~ 15% is necessary for convective carbon 
burning [42]. For lower mass He stars the ^^C abundance was, of course, larger. 
Even with He-star wind mass loss rates reduced by half, Wellstein & Langer 
[10] find a central carbon abundance of > 1/3 at the end of He core burning all 
the way up through 60 Mq stars, so it is clear that convective carbon burning 
will take place. Unfortunately, the cores have not yet been evolved past the CO 
stage. Thus, in the range of ZAMS masses up to > GOMq, if the H envelope 
is lifted off early in the core He burning phase, the convective carbon burning 
will take place after the He burning. 

By ZAMS mass ~ 4OM0, where stars evolve into WR stars almost independent 
of whether they have a companion, the ultimate fate of the compact core is 
uncertain: Brown, Weingartner & Wijers [13] suggest that 1700-37, with a 
progenitor of about 40 Mq went into a low-mass black hole. This would seem 
to indicate that the H-envelope of such massive stars is blown off in an LBV 
phase rapidly enough that the He core again burns as "naked". In any case, 
^^C is burned convectively following He core burning, so the resulting Fe core 
should be small. 

We believe that our discussion earlier in this section indicates that single 
stars in the region of ZAMS masses ~ 20 — 35 Mq end up as high mass 
black holes. We can obtain the high mass black holes, according to our above 
discussion, if we make the He-stars burn with "clothing", i.e., lift their H- 
envelope off only following He core burning. Thus, the evolving massive star 
should meet the companion main sequence star only following He core burning 
(in the supergiant stage). By then its radius R is several hundred Rq, and its 
binding energy 0.6 GM'^/R, very small, because of the large R. In order to 
see effects of matter stripped off from the main sequence companion in the 
transient sources, we want it to end up close to the black hole. Because of its 
low binding energy the supergiant envelope will be expelled by a relatively 
small binding energy of the companion, IMaMbj/c/ where a/ is the distance 
between black hole and companion. In order to make af small, the mass Ma of 
the companion must be small. (More massive main sequence stars will spiral 
in less far, hence end up further from the black hole, and not fill their Roche 
Lobes. However, when they evolve in subgiant or giant phase they will fill it.) 
Both Portegies Zwart, Verbunt, & Ergma [43] and Ergma & Van den Heuvel 
[44] have suggested that roughly the above region of ZAMS masses must be 
responsible for the ~ 7 Mq black holes in the transient X-ray sources in order 
to form enough such sources. Our scenario is essentially the same as that of 
de Kool et al.[45] for the black hole binary A0620— 00. We refer to this work 
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for the properties of the K-star companion, stressing here the evolutionary 
aspects of the massive black hole progenitor. 



6.2 Calculation 



We now calculate the common envelope evolution following the formalism of 
Section 3. Here Ma is the mass of the main sequence companion, M-q that of 
the massive black hole progenitor. The ratio 



is very small and there is great uncertainty in the initial number of binaries 
for such a small q ~ 1/25. We again take the distribution as dq, and again 
assume In a to be uniformly distributed over a logarithmic interval of 7. Again, 
the fraction of binaries in a given interval is 



rfOno^^ (44) 



We evolve as typical a 25Mq star (B) with a companion ~ IMq main sequence 
star (star A) as the progenitor of the transient X-ray sources. The common 
envelope evolution can be done as in Section 3. With M-q^i = 25 Mq and 
neglect of the accretion onto the main sequence mass Ma, we find from Bethe 
& Brown [1] 



i7^"_ 1.2 Mb, (45) 



YiJ ace Ma 

where Y = Mb/ a. Here the coefficient of dynamical friction Cd was taken to 
be 6. The result is relatively insensitive to q, the exponent 1.2 resulting from 
1 + 1/(q-1). 

Thus, in our case 



Yf aceMA\ /0.5Mo\°-^3 , , 

-f = 17 "" ^ =30 — — ^ . 46 

Yi \ Mq I We Ma) ^ ' 

We expect Uce — 0.5, under the assumption that the thermal energy of the 
expelled envelope is equal to that it originally possessed in the massive star 
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(i.e. that it is not used as extra energy to help remove the envelope), but it 
could be smaller. From this we obtain 



a/ M^jYi W,MaJ ' ^ ' 

where we have taken the He star mass Mb,/ to be 1/3 of Mb,!. In order to 
survive spiral-in, the final separation a/ must be sufficient so that the main 
sequence star lies at or inside its Roche Lobe, about 0.2a/ if M\ = Mq. This 
sets ttf ~ 5Rq = 3.5 x lO"*^^ cm and 



/O 5\0-83 
a. = 3.15i^] xlO^^cm, (48) 

\aceJ 

which is about 2 AU. This exceeds the radius of the red giant tip in the more 
numerous lower mass stars in our interval, so the massive star must generally 
be in the supergiant phase when it meets the main sequence star, i.e., the 
massive star must be beyond He core burning. E.g., the red giant tip (before 
the He core burning) for a 2OM0 star is at 0.96 x lO"*^^ cm, for a 25M0 star, 2.5 x 
10^^ cm [46]. These numbers are, however, somewhat uncertain. Notice that 
decreasing ace will increase a,. Decreasing Ma has little inffuence, because with 
the smaller stellar radius the minimum a/ will decrease nearly proportionately. 
Note that neglect of accretion onto the main sequence star would change the 
exponent 0.83 to unity, so accretion is unimportant except in increasing the 
final mass. 

Now a ZAMS 25M0 star ends up at radius 6.7 x lO"*^^ cm (~ 2 Oj) following He 
shell burning [47]. Thus the interval between a^ and 6.7 x 10^^ cm is available 
for spiral-in without merger so that a fraction 



^^^1^0.11 (49) 




3.15 0^ 

of the binaries survive spiral-in, but are close enough so that the main sequence 
star is encountered by the evolving H envelope of the massive star. The He 
core burning will be completed before the supergiant has moved out to ~ 2 
A.U., so binaries which survive spiral-in will have He cores which burn as 
"clothed" , namely as in single stars. 

Given our assumptions in Section 3, the fraction of supernovas which arise 
from ZAMS stars between 20 and 35M0 is 
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where we have assumed the mass lOM© is necessary for a star to go supernova. 
A Salpeter function with index n = 1.5 is assumed here. Our assumption that 
the binary distribution is as dq is arbitrary, and gives us a factor 1/25 for a 
1 Mq companion. Thus, for supernova rate 2 per century, our birth rate for 
transient sources in the Galaxy is 



2 X 10"^ X 0.5 X 0.11 X 0.20 x 0.04 ~ 8.8 x 10"V~^ (51) 

where 0.5 is the assumed binarity, 0.11 comes from eq. (49), and the final (most 
uncertain) factor 0.04 results from a distribution flat in q and an assumed IMq 
companion star. 

In order to estimate the number of transient sources with black holes in the 
Galaxy, we should know the time that a main sequence star of mass ~ ^Mq 
transfers mass to a more massive companion. This depends on the angular- 
momentum loss rate that drives the mass transfer. A guaranteed loss mech- 
anism for close binaries is gravitational radiation, which for a main-sequence 
donor gives a mass transfer rate of lO~^''M0yr~^, almost independent of donor 
mass [48]. As mass is transferred, the mass of the donor decreases and with 
it the radius of the donor. Quite a few low-mass X-ray binaries have X-ray 
luminosities that imply accretion rates in excess of lO~^'^M0yr~^, leading to 
suggestions of additional mechanisms for loss of angular momentum from the 
binary, to increase mass transfer. Verbunt & Zwaan [49] estimate that mag- 
netic braking can boost the transfer of mass in a low-mass binary. We some- 
what arbitrarily adopt an effective mass transfer rate of lO~^M0yr~^ for main 
sequence stars. In order to estimate the number of high-mass black hole, main 
sequence star binaries in the Galaxy we should multiply the birth rate eq. (51) 
by the 10^ yr required, at the assumed mass loss rate, to strip the main se- 
quence star, obtaining 8800 as our estimate. From the observed black-hole 
transient sources Wijers [5] arrives at 3000 low-mass black hole sources in the 
Galaxy, but regards this number as a lower limit. With the uncertainties in 
formation rate and life time, the agreement between the two numbers is as 
good as may be expected. 



6.3 Observations 



We believe that there are many main sequence stars more massive than the 
< 1 Mq we used in our schematic evolution, which end up further away from 
the black hole and will fill their Roche Lobe only during subgiant or giant 
stage. From our earlier discussion, we see that a 2 Mq main sequence star will 
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Table 1 

Parameters of suspected black hole binaries with measured mass functions [5,50- 
54]. N means nova, XN means X-ray nova. Numbers in parenthesis indicate errors 
in the last digits. 



X-ray 


other 
name(s) 


compan. 
type 


(d) 


f{Mx) 

(Me) 


Mopt 

(Mq) 


(/,6) 


names 


q 

(Mopt/Mx) 


Kopt 
(km s-i) 


i 

(degree) 


Mx 

{Mrs) 


d 
(kpc) 


CygX-1 

1956+350 


V1357Cyg 
HDE 226868 


09.7Iab 


5.5996 

74.7(10) 


0.25(1) 


33(9) 
16(5) 


(73.1, +3.1) 

2.5 


LMC X-3 

0538-641 




B3Ve 


1.70 
235(11) 


2.3(3) 


5.6-7.8 


(273.6,-32.1) 

55 


LMCX-1 

0540-697 




07 9III 


4.22 
68(8) 


0.14(5) 




(280.2,-31.5) 

55 


XN Mon 75 

A 0620-003 


V616Mon 
N Mon 1917 


K4 V 
0.057-0.077 


0.3230 

443(4) 


2.83-2.99 

37_44(*) 


0.53-1.22 
9.4-15.9 


(210.0,-6.5) 

0.66-1.45 


XN Oph 77 

H 1705-250 


V2107Oph 


K3 V 


0.5213 

420(30) 


4.44-4.86 
60-80 


0.3-0.6 
5.2-8.6 


(358. 6, +9.1) 

5.5: 


XN Vul 88 

GS 2000+251 


QZ Vul 


K5V 

0.030-0.054 


0.3441 
520(16) 


4.89-5.13 
43-74 


0.17-0.97 
5.8-18.0 


(63.4,-3.1) 

2 


XN Cyg 89 

GS 2023+338 


V404 Cyg 

N Cyg 1938, 1959 


KOIV 
0.055-0.065 


6A7U 

208.5(7) 


6.02-6.12 
52-60 


0.57-0.92 
10.3-14.2 


(73.2,-2.2) 

2.2-3.7 


XN Mus 91 

GS 1124-683 




K5 V 
0.09-0.17 


0.4326 
406(7) 


2.86-3.16 
54-65 


0.41-1.4 

4.6-8.2 


(295.0,-6.1) 

3.0 


XN Per 92 

GRO J0422+32 




MOV 
0.029-0.069 


0.2127(7) 
380.6(65) 


1.15-1.27 

28-45 


0.10-0.97 
3.4-14.0 


(197.3,-11.9) 


XN Sco 94 

GRO J1655-40 




F5-G2 
0.33-0.37 


2.6127(8) 

227(2) 


2.64-2.82 
67-71 


1.8-2.5 
5.5-6.8 


(345.0, +2. 2) 

3.2 


XN 

4U 1543-47 


MX 1543-475 


A2 V 


1.123(8) 
124(4) 


0.20-0.24 
20-40 


1.3-2.6 

2.0-9.7 


(330. 9, + 5. 4) 

9.1(11) 


XN Vol 93 




K6-M0 
0.137± 0.015 


0.2852 
475.4(59) 


3.05-3.29 

~ 78 


0.50-0.65 
3.64-4.74 





(•) A much higher inclination for A0620 has been claimed by Haswell et al. [55] 
of up to i=70. In this case, the lower limits on the component masses would be 
Mx > 3.8 and Mopt > 0.22. 

end up about twice as far from the black hole as the 1 Mq, a 3 Mq star, three 
times as far, etc. Two of the 9 systems in our Table 1 have subgiant donors 
(V404 Cyg and XN Sco). These have the longest periods, 6.5 and 2.6 days 
and XN Sco is suggested to have a relatively massive donor of ~ 2 Mq- It 
seems clear that these donors sat inside their Roche Lobes until they evolved 
off the main sequence, and then poured matter onto the black hole once they 
expanded and filled their Roche Lobe. For a 2 Mq star, the evolutionary 
time is about a percent of the main-sequence time, so the fact that we see 
two subgiants out of nine transient sources means that many more of these 
massive donors are sitting quietly well within their Roche Lobes. Indeed, we 
could estimate from the relative time, that there are 2/9 x 100 = 22 times 
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more of these latter quiet main sequence stars in binaries. 

Amazingly, this factor 22 almost cancels the 1/25 we had for the interval in 
q over which the donors contribute. This is not coincidental. Essentially any 
mass donor, at least almost up to the 25 M© progenitor of the black hole, 
can give rise to a common envelope phase. The BH progenitor crosses the 
Herzsprung gap very quickly, in a time in which the companion can hardly 
accept its mass. (The ratio of g < 1/4 for common envelope evolution was 
determined by Kippenhahn & Meyer-Hofmeister [56] for case A mass transfer.) 
Thus, one can expect essentially all companions, up to g < 1, to go into 
common envelope evolution and contribute. Beginning from Wijers' empirical 
estimate we would thus have (2/9) x 100 x 3000 = 6.7 x 10^ binaries with high- 
mass black holes and main-sequence companions. This number is determined, 
as shown above, chiefly by the number of observed systems with subgiant 
donors. 

If we assume that ZAMS masses ~ 10 — ISM© evolve into a neutron star, we 
should have ~ 3 times more neutron stars than high-mass black holes (see 
eq. (50)). The range follows from our belief that SN 1987A with progenitor 
~ ISMq ZAMS went into a low-mass black hole, following the scenario of 
Brown & Bethe [12]. On the basis of a Monte Carlo calculation using the 
kick velocities of Cordes & Chernoff [21] we find that ~ 1/2 of the binaries 
containing He-star, low-mass main sequence companion (with M ~ IM©) will 
be disrupted in the explosion. Thus, we find only a slightly higher birth rate for 
LMXBs (Low Mass X-ray Binaries) with neutron stars than with black holes, 
although the numbers could be equal to within our accuracy. With comparable 
life times (since the donor masses and mass transfer rates are comparable), 
this would give us one to a few thousand LMXBs with neutron stars, much 
above the total number of observed LMXBs (~ 130). Indeed, from Table 6 of 
Portegies Zwart & Verbunt [57] one sees that their estimated empirical birth 
rate for low-mass X-ray binaries is 2 x 10^^ yr~^, whereas in either theoretical 
evolution including kick velocities they obtain 4 x 10"^ yr^^. This factor of 
20 discrepancy is by far the greatest between theoretical and empirical rates 
in their table, and supports our point that many of the neutron stars must 
have disappeared along the way. Alternatively, a large number of LMXBs with 
neutron stars could be transients as well (like, e.g. AqlX-1). Just at the present 
there are new developments in the evolution of low-mass X-ray binaries, which 
we shall shortly summarize in Section 7. 

As we showed below eq. (48), the He core of the massive star will in general 
be uncovered only after He core burning is completed. The remaining time 
for He burning (in a shell) will be short, e.g., for a 2OM0 ZAMS star it is 
only 1.4 X 10^ years [46]. Therefore the mass loss by wind after uncovering the 
He core will not be large, and when the star finally becomes a supernova, its 
mass will be almost equal to the He core of the original star. The latter can 
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be calculated from 



MHe^O.lO(MzAMs)'-' (52) 

so for ZAMS masses 20 - 35Mq Mhc will lie in the interval ~ 7 - 14 Mq. 

Bailyn et al. [51] find the black hole masses in transient sources to be clustered 
about ~ 7Mq, except for V404 Cyg which has a higher mass. This is in 
general agreement with our scenario, because most of the black holes will 
come from the more numerous stars of ZAMS mass not far from our lower 
limit of ~ 20Mq. Two points are important to note: 

(1) Not much mass can have been lost by wind. Naked He stars have rapid 
wind loss. However in our scenario the He star is made naked only during 
He shell burning and therefore does not have much time (< 10^ yr) to 
lose mass by wind. 

(2) There are good reasons to believe that the initial He core will be rotating 
[58] . The way in which the initial angular momentum affects the accretion 
process has been studied by Mineshige et al. [59] for black hole accretion 
in supernovae. In general accretion discs which are optically thick and 
advection dominated are formed. The disc is hot and the produced energy 
and photons are advected inward rather than being radiated away. The 
disc material accretes into the black hole at a rate of > lO^MEdd for 
the first several tens of days. Angular momentum is advected outwards. 
Our results show that little mass is lost, because the final ~ 7 Mq black 
hole masses are not much less massive than the He core masses of the 
progenitors, and some mass is lost by wind before the core collapses. The 
latter loss will not, however, be great, because there is not much time 
from the removal of the He envelope until the collapse. 

Accretion of the He into the black hole will differ quantitatively from the above, 
but we believe it will be qualitatively similar. The fact that the helium must 
be advected inwards and that little mass is lost as the angular momentum is 
advected outwards is extremely important to establish. This is because angu- 
lar momentum, essentially centrifugal force, has been suggested by Chevalier 
[28] to hold up hypercritical accretion onto neutron stars in common enve- 
lope evolution. (Chevalier [27] had first proposed the hypercritical accretion 
during this evolutionary phase to turn the neutron stars into black holes, the 
work followed up by Brown [3] and Bethe & Brown [1].) However, once matter 
is advected onto a neutron star, temperatures > 1 MeV are reached so that 
neutrinos can carry off the energy. The accreted matter simply adds to the 
neutron star mass, evolving into an equilibrium configuration. Thus, this ac- 
cretion does not differ essentially from that into a black hole. In either case of 
neutron star or black hole an accretion disc or accretion shock, depending on 
amount of angular momentum, but both of radius ~ 10^^ cm, is first formed. 
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giving essentially the same boundary condition for the hypercritical accretion 
for either black hole or neutron star. Thus, the masses of the black holes in 
transient sources argue against substantial inhibition of hypercritical accretion 
by jets, one of the Chevalier suggestions [28]. 

Measured mass functions, which give a lower limit on the black hole mass 
are given in Table 1. Only GRO J0422+32 and 4U 1543-47 have a measured 
mass function < 3Mq. Results of Callanan et al. [60] indicate that the angle 
i between the orbital plane and the plane of the sky for GRO J0422+32 is 
i < 45°, and recent analysis [52] indicate that the angle i for 4U 1543-47 is 
20° <i< 40°. So both GRO J0422+32 and 4U 1543-47 also contain high-mass 
black holes. 

Based on the observations of Kaper et al. [61] that the companion is a hyper- 
giant, Ergma & Van den Heuvel [44] argue that the progenitor of the neutron 
star in 4U1223-62 must have a ZAMS mass > 5OM0. Brown, Weingartner 
& Wijers [13], by similar argumentation, arrived at ~ 45M0, but then had 
the difficulty that 4U1700-37, which they suggested contains a low-mass black 
hole, appeared to evolve from a lower mass star than the neutron star in 1223. 
Wellstein & Langer [10] suggest the alternative that in 1223 the mass occurs 
in the main-sequence phase (Case A mass transfer), which would be expected 
to be quasi conservative. They find that the progenitor of the neutron star 
in 1223 could then come from a mass as low as 26M0. This is in agreement 
with Brown et al. [13] for conservative mass transfer (their Table 1), but these 
authors discarded this possibility, considering only Case B mass transfer in 
which case considerable mass would be lost. 

Wellstein & Langer [10] are in agreement with Brown et al. [13] that 4U1700- 
37 should come from a quite massive progenitor. Conservative evolution here 
is not possible because of the short period of 3.4 days [62]. The compact object 
mass is here 1.8±O.4M0 [63]. Brown et al. [13] suggest that the compact object 
is a low-mass black hole. The upper mass limit for these was found by Brown 
& Bethe [12] to be ~ I.8M0, as compared with an upper limit for neutron star 
masses of ~ 1.5M0. Thus, there seems to be evidence for some ZAMS masses 
of ~ 40 — 5OM0 ending up as low-mass compact objects, whereas we found 
that lower mass stars in the interval from ~ 20 — 35M0 ended up as high-mass 
black holes. In this sense we agree with Ergma & Van den Heuvel [44] that low- 
mass compact object formation "is connected with other stellar parameters 
than the initial stellar mass alone". We suggest, however, following Brown et 
al. [13] that stars in binaries evolve differently from single stars because of 
the different evolution of the He core in binaries resulting from RLOF in their 
evolution. Namely, "naked" He cores evolve to smaller final compact objects 
than "clothed" ones. 

In fact, this different evolution of binaries was found by Timmes, Woosley & 
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Weaver [19]. They pointed out that stars denuded of their hydrogen envelope in 
early RLOF in binaries would explode as Type lb supernovae. They found the 
resulting remnant gravitational mass following explosion to be in the interval 
of 1.2 — 1AMq, whereas in exploding stars of all masses with hydrogen envelope 
(Type II supernova explosion) they found a peak at about I.28M0, chiefly from 
stars of low masses and another peak at 1.73Mq more from massive stars. Our 
Fe core masses in Fig. 1 come from essentially the same calculations, but the 
"Remnant" masses of Woosley & Weaver [18] are somewhat greater than those 
used by Timmes et al. [19]. In fact, the differences between the masses we plot 
and those of Timmes et al. come in the region ~ 1.7 — I.SMq (gravitational). 
This is just in the Brown & Bethe [12] range for low-mass black holes. It 
may be that some of the stars with low-mass companions evolve into low- 
mass black holes. Presumably these would give lower luminosities than the 
high-mass black holes, although at upper end of the mass range we discuss 
4U1700-37 seems to be an example of such a system. Of course here the high 
luminosity results from the high mass loss rate of the giant companion. There 
are substantial ambiguities in fallback, etc., from the explosion. Our point in 
this paper is that most of the higher mass single stars 20 — 35Mq go into 
high mass black holes. (The Brown & Bethe [12] limit for low-mass black hole 
formation is ~ 1.5 — I.SMq gravitational, but there is some give and take in 
both lower and upper limit. Also the stars are not all the same. In particular 
different metallicities will give different wind losses.) 



7 Evolution of Low Mass X-ray Binaries 



We shall briefly point out new developments in the evolution of low-mass X- 
ray binaries. These were foreseen in the excellent review by Van den Heuvel 
[64], and there has been substantial development in this field lately. 

Low mass X-ray binaries are considered to be progenitors of recycled pulsars 
with helium white dwarf companions. In order to bring the magnetic fields 
of the latter down to ~ 10^ G and to speed them up to their final period. 
Van den Heuvel & Bitzaraki [65] had the neutron star accreting ~ 0.5 Mq 
from the main-sequence donor. More detailed recent calculations by Tauris & 
Savonije [66] find that if the initial orbital period is below ~ 30 days with a 
main sequence donor of ~ 1 Mq which undergoes stable mass transfer with 
the neutron star, the mass of the latter is increased up to ~ 2 Mq if the 
amount of material ejected as a result of propeller effect or disk instabilities 
is insignificant. This presents a problem for us because the Brown & Bethe 
[12] mass limit for neutron stars is 1.5 Mq. From this limit, we would say that 
these neutron stars in low-mass X-ray binaries would have gone into black 
holes. 
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A way out of this problem was suggested by Van den Heuvel [64], which is 
called the evolution of Her X-1 type X-ray binaries (see especially the Ap- 
pendix of Van den Heuvel [64]). In this case a radiative donor more massive 
than the neutron star pours matter in unstable mass transfer across the Roche 
Lobe onto the neutron star. This mass transfer can occur onto the accretion 
disc by as much as ~ 10"^ M^dd, if ^Edd ~ 1-5 x 10~^ Mq yr~^ is accreted 
onto the neutron star, since the Eddington limit goes linearly with R and the 
radius of the disc can be ~ 10^'^ cm. The advection dominated inflow-outflow 
solution (ADIOS) of Blandford & Begelman [23] suggests that the binding en- 
ergy of the matter released at the neutron star can carry away mass, angular 
momentum and energy from the gas accreting onto the edge of the accretion 
disc provided the latter does not cool too much. In this way the binding energy 
of a gram of gas at the neutron star can carry off ~ 10^ grams of gas at the 
edge of the accretion disc. Such radiatively-driven outflows are suggested by 
King & Begelman [22] to enable common envelope evolution to be avoided. 
Tauris & Savonije [66] have carried out a detailed evolution of low-mass X-ray 
binaries with Porb > 2 days using computer programs based on Eggleton's, 
which for radiative and semiconvective donors follow, in at least a general 
way, the above ideas. For a deeply convective donor a short phase of rapid 
mass loss may reach a rate as large as lO^MEdd while the mass of the donor 
drops to well below the neutron star mass. Although rates > 10^ Mgdd would 
be hypercritical for spherical accretion, somewhat higher rates survive hyper- 
critical accretion provided angular momentum is taken into account [28] . The 
important point is that the donor mass can be brought down sufficiently far 
before stable mass transfer at a rate < M-^dd sets in, so that the neutron star 
can avoid accreting sufficient mass to send it into a black hole. It is not clear 
what percentage of the neutron stars will survive black-hole fate. Our rough 
estimates in Section 6 indicate that only a small fraction need to do so. 

For even more massive donors (2-6 Mq) which are either radiative or semicon- 
vective, work by Tauris et al. [24] indicates that the low-mass X-ray binaries 
with C/0 white dwarf (CO) companions can be made in much the same way. 
In an earlier paper. Van den Heuvel [67] had suggested that these binaries 
would originate from donor stars on the asymptotic giant branch. In order to 
evolve these, he needed an efficiency a > 1; i.e., sources additional to those 
included in our earlier common envelope evolution, such as mass loss by in- 
stabilities in the AGB, dissociation energy, etc. have to participate in helping 
to remove the envelope of the donor star. 

King and Ritter [25] have computed a scenario for CygX-2 with an initial 
donor mass of ~ 3.6Mq. Currently the donor has a mass of 0.5 — 0.7Mq and a 
large radius, about 7Rq. About 2Mq must have been lost in super-Eddington 
accretion, roughly along the lines sketched above. More massive donors can 
lead to relatively more massive white dwarf companions, which will be C/0 
white dwarfs. 
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In fact, the present situation is that no circular NS-COcQ binaries which 
went through common envelope evolution seem to be observed, the alterna- 
tive Tauris et al. [24] evolution which avoids common envelope evolution being 
preferred. This presents a real dilema for the standard scenario of common 
envelope evolution. It seems clear [68] that in the binary B2303+46 the com- 
panion to the pulsar is a C/0 white dwarf. B2303+46 is an eccentric binary 
NS-COe, indicating that the neutron star was formed last. This is confirmed 
by the unrecycled field strength of the pulsar B = 8 x 10^^ G. Cases have 
made that the recently discovered J1141— 6545 [69] and B1820— 11 [70] are 
also NS-COe binaries. 

On the other hand, evolutionary calculations show that formation probability 
of NS-COc binaries through common envelope evolution is > 50% as probable 
as of NS-COe binaries [71]. In this evolution the pulsar magnetic moment 
will be recycled, brought down at least a factor of 100 [3] and possibly even 
further, down to the empirical values of ~ 5 x 10*^ G found in the NS-COc 
binaries. The lowering of the magnetic fields increases the time of observation 
by a factor of ~ 100 or of 2000, depending on whether the theoretical or 
empirical magnetic field is used. Since we fairly certainly observe at least one 
NS-COe binary, we should see either 100 or 2000 NS-COc binaries which have 
gone through common envelope evolution. We certainly don't see anything 
like this, at most the 5 that had earlier been attributed to common envelope 
evolution, and probably none. Brown, Lee, & Portegies Zwart [71] remove 
at least most of this discrepancy by showing that with the introduction of 
hypercritical accretion the neutron star in common envelope evolution with 
the evolving main sequence companion goes into a black hole. 



8 Discussion and Conclusion 



Our chief new point in the evolution of binaries of compact objects is the use 
of hypercritical accretion in common envelope evolution, although the idea 
of hypercritical accretion is not new (Section 3). Chevalier [28] discussed the 
possibility that angular momentum might hinder hypercritical accretion. In 
his treatment of the accretion disc, he assumed gas pressure to dominate, in 
order to raise the temperature sufficiently for neutrinos to be emitted. This 
entailed a tiny viscosity, characterized by a < 10~^ in the a-description. More 
reasonable values of a are ~ 0.1. 

Bethe, Brown, & Lee [72] have shown that for larger a's, a ~ 0.01 — 1, the disc 
pressure is radiation dominated, and they find a simple hypercritical advection 
dominated accretion fiow (HADAF) of matter onto the neutron star. 



The lower suffix c (e) denotes the circular (eccentric) binaries. 
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The Bethe, Brown, & Lee HADAF appears to reproduce the Armitage & Livio 
[77] numerical two-dimensional hydro solution. These latter authors suggest 
that jets will prevent hypercritical accretion by blowing off the accreting mat- 
ter. At such high rates of accretion ~ IMq yr~^ the Alfven radius is, however, 
close to the neutron star surface, and we believe that this will effectively, shut 
down any magnetically driven jets. 

In Section 7 we discussed the advection of a rotating He envelope into a black 
hole. We believe that two possibilities exist. Phinney & Spruit suggest [80] 
that the magnetic turbulence is strong enough to keep the He envelope in 
corotation with the core of the star until shortly before it evolves into a black 
hole. Then not much angular momentum would have to be advected away in 
order to let the matter accrete. Alternatively, magnetic turbulence is strong 
enough so that angular momentum can be carried away from a rapidly rotating 
He core; then the matter can accrete. From the measured masses of ~ 7 Mq 
we know that most of the He core must fall into the black hole, so one of these 
scenarios should hold. Both favor high magnetic turbulence, lending credence 
to the Chevalier suggestion we quoted. 



Appendix 

A Common Envelope Evolution of Cygnus X-3 



The closeness of the compact object in CygX-3 to its ~ 10 Mq compan- 
ion helium star bears witness to an earlier stage of common envelope evo- 
lution. Although the mass of the He star has not been measured, the star 
is similar to V444 Cygni, the mass of which is 9.3 ± 0.5 Mq [81]. For ex- 
ample, from the period change its mass loss rate would be M^y^ = 0.6 x 
10-5(MHe/10 Mq) Mq yr-i [82] whereas that of V444 Cygni is Mdyn = 
1 X 10~^ Mq jt"^ [83] indicating an Mhc ~ 10 Mq. Mass loss rates cannot 
easily be obtained from W.-R. winds because of large nonlinear effects which 
necessitate corrections for "dumpiness". However polarization measurement 
of the Thomson scattering, which depend linearly on the wind, give a mass 
loss rate of ~ M = 0.75 x 10~^ Mq yr~^ [84], roughly compatible with the 
period change. In agreement with many other authors we take M^e = 10 Mq 
in CygX-3. 

Here we evolve a massive 0-star binary with initial ZAMS masses of 33 Mq 
and 23 Mq as possible progenitor for CygX-3. In red giant phase the 33 Mq 
star will transfer its H envelope to the 23 Mq companion, leaving a He star of 



35 



MHe = 0.1 M^AMS = 13 Mq. (A.l) 

With efficiency of mass transfer assumed to go as g^, about half of the 20 Mq 
H-envelope will be accepted by the companion, which then becomes a reju- 
venated 33 Mq star. The He core of the primary then explodes, going into a 
1.5 Mq compact object, neutron star or low-mass black hole. After the com- 
panion 33 Mq star evolves, the binary will go into common envelope evolution. 
Eq. (29) can be written 






where we again take q = 6. With Mb,^ = 33 Mq and MA,i = 1.5 Mq, 

Yf/Y, = 27. (A.3) 

The compact object mass scales as 

Ma oc rV(cd-i) = yi/5 ^A^4) 

so that 

Ma,/ = 2.9 Mq (A.5) 

and the final compact object is certainly a black hole, in agreement with 
Cherepaschchuk & Moffat [85] and Ergma & Yungelson[86]. We believe our 
evolution here to show that this ~ 3 Mq black hole is about the most massive 
that can be formed in common envelope evolution by accretion onto a low- 
mass compact object, since our 33 Mq companion is near to the ZAMS mass 
range that will lose mass in an LBV phase, unsuitable for common envelope 
evolution, so it cannot be made much more massive. We next find 

tti Mn i Yf . , X 

-1. = — Mii ^ ~ 70. A.6 

a/ Mb,/ Y, ^ ' 

For an a/ ~ 3.5 Rq this gives 

tti ~ 250 Rq (A.7) 

comfortably within the red giant range. 
Following Ergma & Yungelson [86] we calculate the accretion rate as 
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M„o = 0.14 (^)'«r„to P^" (^)'"^fw,„a. (A.8) 

Here f looo is the wind velocity in units of 1000 km s~^ and Phr is the period 
in hours.0 For Mwind we, as Ergma & Yungelson, take M^yn- These authors 
take fiooo = 1-5, essentially the result of Van Kerkwijk et al. [82]. An earlier 
estimate by Van Kerkwijk et al. [87] was Viooo = 1- We believe that the fwind 
to be used here may be different from the (uncertain) measured terminal wind 
velocities, because the velocity near the compact object is substantially less. 
Therefore, we take t^iooo = 1- Taking Mwind = Miyn we obtain 



Mace = 2.2 X 10-' Mo yr-i (A.9) 

This is to be compared with 

MEdd = ^ncR/n,, = 2.6 x 10^^ (Mbh/M©) Mq yr'^ (A. 10) 

where k^s = 0.2 g/cm for He accretion. Our result is in fair agreement with 
Ergma & Yungelson [86], who find Msdd ~ 2.3 x 10"' Mq yr"^ for a 10 M© 
black hole. The presence of jets in CygX-3 argues for super Eddington rates 
of accretion, which we find. 

Cherepaschchuk and Moffat [85] estimated the total luminosity of CygX-3 to 
be Lboi ~ 3 X 10^^ erg. The efficiency of black hole accretion varies as 

0.057<e<0.42 (A.ll) 

for a black hole at rest to a (maximally rotating) Kerr black hole. We expect 
the black hole to be spun up by accretion from the wind or accretion disc. 
Taking an intermediate e = 0.2, we find 

L = 2.5 X 10^^ erg s"^ (A.12) 

in rough agreement with the Cherepaschchuk and Moffat value. 

CygX-3 is often discussed as the "missing link" in binary pulsar formation. 
In fact, because of its high He star mass, upon explosion of the latter, it most 
probably will break up. But it should be viewed as "tip of the iceberg" [64], in 
that there must be a great many more such objects with lower mass He stars 



^ Through a slip, the two factors preceding M^ind appear in the denominator in 
[86], although we confirm that they carried out their calculations with the correct 
formula. 



37 



which are not seen. We have shown, in section 3 however, that these objects 
are more hkely to contain a black hole than a neutron star. 

In our evolutionary scenario, the He star progenitor has about the same ZAMS 
mass as that of the primary. Thus, the fate of the "naked" He star should be 
the same low-mass compact object, neutron star or low-mass black hole that 
resulted from the explosion of the primary. 



B Implications for LIGO 



Our results that there are 10 times more [] black hole, neutron star binaries 
than binary neutron stars has important results for LIGO, the detection rates 
of which were based on the ~ 10~^ per year per galaxy rates of merging for 
the latter. The combination of masses which will be well determined by LIGO 
is the chirp mass 



Mchirp = ^?'^M^"' = {MiM^f"' (Ml + Ma)-^/' (B.l) 

where M = Mi + M2 is the total system mass. The chirp mass of a NS-NS 
binary, with both neutron stars of mass 1.4 Mq, is 1.2 Mq. A 10~^ birth rate 
implies a rate of 3 yr~^ out to 200 Mpc [88]. Kip Thorne informs us that 
LIGO's first long gravitational- wave search in 2002-2003 is expected to see 
binaries with Mchirp = 1-2 Mq out to 21 Mpc. 

The chirp mass corresponding to the Bethe & Brown [1] LMBH-NS binary 
with masses 2.4 Mq and 1.4 Mq is 1.6 Mq. Including an ~ 30 % increase in 
the rate to allow for high-mass black-hole, neutron-star mergers (which should 
be regarded as a lower limit because of the high-mass limit of 80 Mq used by 
Bethe & Brown for going into a HMBH) gives a 26 times higher rate than 
Phinney's estimate for NS-NS mergers. There factors are calculated from the 
signal to noise ratio, which goes as MJ^-^^.^ and then cubing it to obtain the 
volume of detectability. We then predict a ratio of 3 x (21/200)^ x 26 = 0.09 
-yj,-i £qj, 2003, rather slim. The enhanced LIGO interferometer planned to 
begin in 2004 should reach out beyond 150 Mpc for Mchirp = 1-2 Mq, increasing 
the detection rate to 3 x (150/200)^ x 26 = 33 yr~^. We therefore predict that 
LIGO will see more mergers per month than NS-NS mergers per year. 



^ Actually about 20 times more if we include the binaries in which the pulsar goes 
into a black hole in the He shell burning evolution. However, these will have masses 
not very different from the binary neutron stars so we do not differentiate them. 
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C Binary Contributions to Gamma Ray Bursters 



The sheer numbers of black-hole, neutron-star binaries should dominate the 
mergers for gravitational waves, which could be detected by LIGO [39]. For 
gamma-ray bursts, the presence of an event horizon eases the baryon pollution 
problem, because energy can be stored in the rotational energy of the black 
hole, and then released into a cleaner environment via the Blandford-Znajek 
magnetohydrodynamic process. 

Binaries containing a black hole, or single black holes, have been suggested 
for some time as good progenitors for gamma-ray bursts [89-92,34,93]. Rea- 
sons for this include the fact that the rest mass of a stellar mass black hole 
is comparable to what is required to energize the strongest GRB. Also, the 
horizon of a black hole provides a way of quickly removing most of the mate- 
rial present in the cataclysmic event that formed it. This may be important 
because of the baryon pollution problem: we need the ejecta that give rise to 
the GRB to be accelerated to a Lorentz factor of 100 or more, whereas the 
natural energy scale for any particle near a black hole is less than its mass. 
Consequently, we have a distillation problem of taking all the energy released 
and putting it into a small fraction of the total mass. The use of a Poynting 
flux from a black hole in a magnetic field [94] does not require the presence 
of much mass, and uses the rotation energy of the black hole, so it provides 
naturally clean power. 

As a neutron star in a binary moves nearer to a black hole companion, it is 
distorted into a torus around the latter. Most of the torus matter enters the 
black hole from the last stable Keplerian orbit oi R = 6GMbh/c^, carrying 
considerable angular momentum. In the process the black hole is spun up 
until it rotates with some fraction of the speed of light. A magnetic field 
which originates from the neutron star, but which could have been enhanced 
by differential rotation is anchored in the remaining part of the torus, the 
accretion disc. 

When a rapidly rotating black hole is immersed in a magnetic field, frame 
dragging twists the field lines near the hole, which causes a Poynting flux to 
be emitted from near the black hole. This is the Blandford-Znajek mechanism 
[94] . The source of energy for the flux is the rotation of the black hole. The 
source of the field is the surrounding accretion disk or debris torus. We showed 
[95] that at most 9% of the rest mass of a rotating black hole can be converted 
to a Poynting flux, making the available energy for powering a GRB 

Ebz = 1.6 X 10^^ (M/Mq) erg. (C.l) 

The power depends on the applied magnetic field: 
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Pbz - 6.7 X 10^° 5^5 {M/MqY erg s"^ (C.2) 

(where Bi^ = B/10^^ G). This shows that modest variations in the apphed 
magnetic field may explain a wide range of GRB powers, and therefore of 
GRB durations. There has been some recent dispute in the literature whether 
this mechanism can indeed be efficient [96] and whether the power of the BH 
is ever significant relative to that from the disk [97]. The answer in both cases 
is yes, as discussed by Lee, Wijers, & Brown 



The issue, therefore, in finding efficient GRB sources among black holes is 
to find those that spin rapidly. There are a variety of reasons why a black 
hole might have high angular momentum. It may have formed from a rapidly 
rotating star, so the angular momentum was there all along ('original spin', 
according to Blandford [98]); it may also have accreted angular momentum 
by interaction with a disk ('venial spin') or have formed by coalescence of a 
compact binary ('mortal spin'). We shall review some of the specific situations 
that have been proposed in turn. 

Neutron star mergers are among the oldest proposed cosmological GRB sources 
[99-101], and especially the neutrino fiux is still actively studied as a GRB 
power source [102]. However, once the central mass has collapsed to a black 
hole it becomes a good source for BZ power, since it naturally spins rapidly 
due to inheritance of angular momentum from the binary [103]. Likewise BH- 
NS binaries [104] will rapidly transfer a large amount of mass once the NS fills 
its Roche lobe, giving a rapidly rotating BH [105]. The NS remnant may then 
be tidally destroyed, leading to a compact torus around the BH. It is unlikely 
that this would be long-lived enough to produce the longer GRB, but perhaps 
the short (t < 1 s) ones could be produced [106]. However, mass transfer could 
stabilize and lead to a widening binary in which the NS lives until its mass 
drops to the minimum mass of about 0.1 Mq, and then becomes a debris torus 
[107] By then, it is far enough away that the resulting disk life time exceeds 
1000 s, allowing even the longer GRB to be made. Thus BH-NS and NS-NS 
binaries are quite promising. They have the added advantage that their envi- 
ronment is naturally reasonably clean, since there is no stellar envelope, and 
much of the initially present baryonic material vanishes into the horizon. 

In addition to the mergers from compact objects. Fryer & Woosley [34] sug- 
gested that GRBs could originate from the coalescence of low-mass black hole 
and helium-star binaries in the Bethe & Brown [1] scenario. From eq. (35) 
we see that binaries survived in the initial range of 0.5 x lO^^cm < a^ < 
1.9 X lO^^cm. Inside that range for 0.04 x lO^^cm < a^ < 0.5 x 10^^ cm the 
low-mass black hole coalesces with the core. Hence, using a separation distri- 
bution fiat in In a, coalescences are more common than low-mass black-hole, 
neutron-star binaries by a factor ln(0.5/0.04)/ln(1.9/0.5) = 1.9. In Bethe & 
Brown [1] the He star compact-object binary was disrupted ~ 50 % of the 
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Table 2 

Summary of the formation rates of various sources of gamma-ray bursts (GRB) or 
gravity waves (GW) from the binaries considered in this review. L(H)BH means 
low- (high-)mass black hole. 

object GRB GW rate [GEM*] 



NS - NS merger XX 10 

NS - BH merger X X 100 

WR star - LBH merger X 380 



* GEM means Galactic Events per Megayear; rates are quoted for redshift 0. 

time in the last explosion, which we do not have here. Thus, the rate of low- 
mass black-hole. He-star mergers is 3.8 times the formation rate of low-mass 
black-hole, neutron-star binaries which merge, or R = 3.8 x 10~^ yr~^ in the 
Galaxy. 

In Table 2 we summarize the formation rates of GRBs and gravity waves from 
the binaries considered in this review. 

Because gamma-ray bursts have a median redshift of 1.5-2 [108]), and the su- 
pernova rate at that redshift was 10-20 times higher than now, the gamma-ray 
burst rate as observed is higher than one expects using the above rates. How- 
ever, for ease of comparison with evolutionary scenarios we shall use the GRB 
rate at the present time (redshift 0) of about 0.1 GEM. (Wijers et al. [108] 
found a factor 3 lower rate, but had slightly underestimated it because they 
overestimated the mean GRB redshift; see ref. [106] for more extensive discus- 
sions of the redshift dependence). An important uncertainty is the beaming 
of gamma-ray bursts: the gamma rays may only be emitted in narrow cones 
around the spin axis of the black hole, and therefore most GRBs may not be 
seen by us. An upper limit to the ratio of undetected to detected GRB is 600 
[109], so an upper limit to the total required formation rate would be 60 GEM. 
We may have seen beaming of about that factor or a bit less in GRB 990123 
[110], but other bursts (e.g. 970228, 970508) show no evidence of beaming 
in the afterglows (which may not exclude beaming of their gamma rays). At 
present, therefore, any progenitor with a formation rate of 10 GEM or more 
should be considered consistent with the observed GRB rate. 

An exciting possibility for the future will be to receive both gravitational- wave 
and gamma-ray burst signals from the same merger, with attendant detailed 
measurement, which would give witness to them arising from the same binary. 

Because we dealt in this review with binaries, we did not explain one popular 
model of GRBs, the Woosley Collapsar model [92] . In this model a black hole 
is formed in the center of a rotating W.-R. star. The outer matter can then be 
accreted into the neutron star, spinning it up. If, however, magnetic turbulence 
is sufficient to keep the envelope of the progenitor in corotation with the core 
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until a few days before collapse of the latter, as suggested by Phinney & Spruit 
[80] the He envelope could not furnish enough angular momentum to the black 
hole for the latter to drive the necessary jets (see end of Section 8). 
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